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Abstract — In this paper, we consider the downlink signal-to- 
interference-plus-noise ratio (SINR) analysis in a heterogeneous 
cellular network with K tiers. Each tier is characterized by 
a base-station (BS) arrangement according to a homogeneous 
Poisson point process with certain BS density, transmission 
power, random shadow fading factors with arbitrary distribution, 
arbitrary path-loss exponent and a certain bias towards admitting 
the mobile-station (MS). The MS associates with the BS that 
has the maximum SINR under the open access cell association 
scheme. For such a general setting, we provide an analytical 
characterization of the coverage probability at the MS. 

Index Terms — Multi-tier networks, Cellular Radio, Co-channel 
Interference, Fading channels, Poisson point process. 

I. Introduction 

THE heterogeneous cellular network is a complex overlay 
of multiple cellular communication networks such as the 
macrocells, microcells, picocells, and femtocells. Research has 
shown heterogeneous networks support greater end-user data- 
rate and throughput as well as better indoor and cell-edge 
coverage. This has further led to its inclusion as an important 
feature for the 4G cellular networks |Q]-E]. 

In the conception of the heterogeneous cellular network, 
one is looking at an overlay of several dense, irregularly 
and often completely randomly deployed networks (namely 
microcells, picocells and femtocells) with a limited cover- 
age area, all deployed on top of the conventional macrocell 
network. These networks consist of base-stations (BSs) with 
different transmission powers, different traffic-load carrying 
capability, and different radio environment which is based on 
the locations in which they are deployed. All these sum up to 
an extremely complicated network. As a result, the analysis of 
such a network through system simulations (which is largely 
the approach taken for studying the conventional macrocell 
network) is hampered by the curse of dimensionality due to 
the many parameters involved in designing and modeling each 
of the representative networks that make up the heterogeneous 
network. For this reason, we seek to develop an analytical 
model that captures all the design scenarios of interest. 

The analysis in this paper applies to the downlink perfor- 
mance in terms of the signal-to-interference-plus-noise ratio 
(SINR) at the mobile-station (MS), where the MS associates 
itself with the BS that has the maximum SINR at the MS. 
The SINR at the MS is an important metric that determines 
the outage probability (or coverage probability), capacity and 
throughput of a cellular network in the downlink, and the char- 
acterization of the distribution of SINR aids in the complete 
understanding of the SINR metric. 
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The recent focus on stochastic geometry as a means to 
study and analyze large systems of essentially randomly 
deployed nodes, combined with the fact that certain rep- 
resentative networks (femtocell networks), that are a part 
of the heterogeneous network, are formed due to the end- 
user deployments, thus falling under the random deployment 
category, has naturally led to the application of stochastic 
geometry in modeling the heterogeneous network. While there 
is a vast literature corresponding to stochastic geometric 
modeling and analysis of systems with randomly deployed 
nodes, there have been three independent efforts to applying 
stochastic geometry to study heterogeneous networks ||5l- 
171 . All the 3 papers modeled the heterogeneous network 
as being composed of multiple tiers where the BSs in each 
tier are deployed according to a homogeneous Poisson point 
process, independent of the other tiers. While the first two 
references focused on the case where the fading coefficients 
were modeled as independent and identically distributed (i.i.d.) 
random variables with exponential distribution and obtained 
closed-form expressions for the distribution of SINR at the 
MS in a heterogeneous network for positive values of SINR 
(in dB), our work has provided semi-analytical expressions for 
the same quantities and holds for arbitrary fading distributions 
and for all values of SINR. 

Jo et. al. [8 1 have extended their results to obtain a complete 
characterization of the SINR in the heterogeneous network 
where the exponents for the power-law path-loss model were 
different for different tiers of the heterogeneous network, but 
for the case when the MS associates itself with the nearest 
BS rather than the BS that maximizes the SINR, and when 
an exponential distribution was assumed for fading. With this 
they could study the effect of varying cell-association biasing. 

In this paper, we incorporate all the modeling details avail- 
able in the literature for studying the heterogeneous cellular 
network, introduce a few additional important features, and 
obtain the complete characterization of the downlink SINR, 
and the coverage probability (i.e. 1 - outage probability) in 
a heterogeneous network, for the original case where the MS 
associates with the BS with the maximum SINR at the MS. 
In particular, the following paragraph lists our contributions. 

From the modeling stand-point, the transmission and chan- 
nel characteristics of a BS corresponding to a given tier 
include the transmission power, bias factor, path-loss exponent, 
distribution of fading coefficients, and the SINR threshold 
to be satisfied by the MS in order to communicate with 
the BS. These may be different for different tiers. This is 
more general than the setting in the previous heterogeneous 
network analysis literature. For this situation, we have obtained 
accurate characterization for the SINR distribution as well as 
the coverage probability for a given MS. Since we are able to 
handle arbitrary fading distributions, that are further different 
for different tiers, we are able to consider practical fading 
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models and obtain the coverage probabilities in those cases. 
Such a strong result is likely to be useful for studying and 
analyzing realistic scenarios. The following section introduces 
the system model. 

II. System Model 

This section describes the various elements used to model 
the wireless network, namely, the BS layout, the radio envi- 
ronment, and the performance metrics of interest. 

1) BS Layout: The BS layout for the k th tier, where 
k G {1, . . . , K} , is according to independent homogeneous 
Poisson point process with density, A& . 

2) Cell-Association policy: The MS associates itself to the 
BS corresponding to the the strongest instantaneous received 
power (the BS from which the MS has the maximum SINK). 
We focus on the open access scheme in this paper where the 
MS can freely communicate with any of the K tiers. 

3) Radio Environment: The received power at the MS 
from the j th BS belonging to the k th tier at a distance 
D kj (> 0) from the MS is given by P = I h, ■ 
where (Pk, ^kj, Bk, £fc) corresponds to the constant trans- 
mission power, random channel gain coefficient, constant bias 
coefficient and the constant path-loss exponent (> 2) of the 
k th tier, respectively. Further, ^kj can assume any arbitrary 
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< oo. It is independent and 



distribution as long as E 

identically distributed (i.i.d.) across all the BSs of the fc th tier, 
and independent of the other tiers and the underlying random 
process governing the BS arrangement. For this reason, we 
drop the subscript j and for compactness denote the expecta- 
tion by . In essence, the different tiers have their own 
arbitrary fading distributions. All BSs of the k th tier adopt 
an identical bias factor, Bk, which is some positive value. By 
manipulating the Bk's corresponding to the different tiers, we 
may regulate the traffic from one tier to the other. 

4) Performance Metric: In this paper, we are concerned 
with the SINR at a given MS. Without loss of generality, the 
MS is assumed to be located at the origin of the plane. SINR 
is defined as the ratio of the received signal from the desired 
BS to the sum of the interferences from all the BSs belonging 
to all the K tiers and the background noise. As a result 



SINR = 



Bk 



Em=l X^=l^ro*7 
(m,l)^(k, j) 



i.RJ] m B r , 



(1) 



V 



where (k, j) corresponds to the indices of the tier and the 
corresponding BS, which has the maximum received power at 
the MS, and rj is the power corresponding to the background 
noise. Further, each user can successfully communicate with 
its desired BS provided the SINR is above the minimum 
threshold value, /3k, that is a characteristic of the tier. As a 
result, the coverage probability is defined as the probability 
that the MS is able to communicate with a BS using a specific 
BS association protocol. 

III. Useful Lemmas 

We will first present some lemmas which will be useful in 
deriving the coverage probability. 



Lemma 1. The SINR at the MS has the same distribu- 
tion as that of a single-tier network where all BSs in the 
network have unity transmission power, channel gain and 
path-loss exponent, and are arranged according to a non- 
homogeneous 1-D Poisson point process with BS density 

i , r > 0, as 



function A(r) = YaLi x i (PiBi)^ E*, S! rn \ 

_2_ 

long as E*; 1 < oo, V I = 1, 2, ■ • • , K. That is 



SINR 
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A(r) 



where = st indicates the equivalence in distribution; and 
j is the ascendingly ordered distances of the BSs from 
the origin, obtained from a non-homogeneous 1-D Poisson 
point process with BS density function A (r) defined above. 

Proof: See Appendix lAl ■ 
Now, for the equivalent single-tier network, we characterize 
the distance of the nearest k th tier BS from the MS, i^" 

Lemma 2. The tail probability of R^' is 

P({^ fe) >r}) 
= exp ( -A fc 7r(P feJ Bfc)^ E *A 



, V r > 0. (3) 



Proof: See Appendix IB"1 ■ 
Using the above result, we will now characterize the random 
variable /, which represents the tier to which the desired BS 
(i.e., the BS nearest to the MS in the equivalent single-tier 
network mentioned in Lemma [TJ belongs. 

Lemma 3. The desired BS belongs to the fc th tier 



(fc = l, 2,- 



K) with the probability 
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t~ dt. (4) 



Further, in the special case {efcj-^j = £ 



'({I = k}) = 



A fe E*= {P k B k )' 
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(5) 



Proof: See Appendix ICl ■ 
Although we do not have a closed-form expression, (|4]) 
can be computed easily to any desired accuracy by numerical 
integration. Next, we obtain the p.d.f. of the distance of the 
serving BS from MS, given that it belongs to the k th tier, 
denoted by f^\j (r\k). 

Lemma 4. The p.d.f. of the distance of the serving BS from 
the MS, given it belongs to the fc th tier (fc = 1, 2, • • • , K) , 
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P ({/ = *}) 

xexp ^-5^A I 7r(flB,)^E#f r r^J , r > 0, (6) 

where the denominator is obtained using Lemma [3] 

Proof: See Appendix |B"1 ■ 
Now, we are ready to obtain the expression for the coverage 
probability for a typical MS in this heterogeneous network. 

IV. Coverage Probability 

Recall that the MS is covered only if the SINR exceeds a 
certain threshold, f3 k , where fc is the tier to which the desired 
BS belongs. Using the results from the previous section, 
we now present the coverage probability for an MS in a 
heterogeneous cellular network. 

Theorem 1. The coverage probability of the MS is 
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Proof: See Appendix [E] ■ 
Though the above expression is not in closed form, this 
can also be computed to any desired accuracy using numerical 
integration methods to compute the double integral. Moreover, 
several insightful results arise for certain special cases of the 
above result. For example, in the interference-limited scenario, 
where the effect of the background noise may be ignored in the 
presence of the strong interferences from the BSs, the above 
result reduces to a simple form as shown below. 

Corollary 1. In the interference-limited case, with {s k }^ =1 = 
e, the coverage probability is 



rniopen— access 
coverage 



A fe E*= (P k B k ) 



Ik 



K 
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where lk = f™ = _ 



duj, and this integra- 



tion reduces to 7 fc = sm ^j/ £) • /3 k e , if Pk > E 

The detailed proof will not be given. However, the corollary 
can be easily proved by noting that the integration w.r.t. 

t in (1711 evaluates to — -, — 5 5 — r-^- — . Notice that 

7^ does not depend on any of the parameters that define 
the characteristics of the BSs of the various tiers, and only 
depends on the path-loss exponent, and the SINR thresholds 
of the various tiers. Further, j k is the coverage probability 
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Figure 1. 2-Tier network, Coverage Probability vs Tier 1 SINR threshold 

of a MS in a single-tier network. For the 77 = case, J5] 
Eq. 3] gives the coverage probability of a single-tier network 
with SINR threshold, 7 (> 1) , and channel gains that are 
i.i.d. exponentially distributed to be s '" 2 ^^ £ ' ) • 7~?. For the 
same case, in [9 Remark 4], we have shown that the single- 
tier network SINR distribution is the same irrespective of 
the transmission power of BSs and the distribution of the 
channel gains. As a result, the above expression that holds 
for exponential fading distribution also holds for any other 
fading distribution (even no fading). 

The contributions of the BS density, transmission power, 
fading distribution, and the bias factor of the k th tier BS, 
are all captured as the multiplicative factors of the k th tier 
SINR threshold, j k . Further, the contribution of the fading 
coefficient is completely captured in terms of the ^- moment 
of the random variable ^ k , used to represent the shadow 
fading factor of the k th tier. Next, having studied the coverage 
probability of the MS in the heterogeneous cellular network, 
we present some numerical examples to illustrate them. 



V. Numerical Examples and Discussion 

In this section, we study various scenarios in order to clearly 
illustrate the results we have obtained. We restrict ourselves 
to a two-tier network consisting of a macrocell and picocell 
network for simplicity, and assume that the background noise 
is zero. We note that these studies can be extended to arbitrary 
number of tiers. Further, please refer Appendix for the 
algorithm used to perform the Monte-Carlo simulations. For all 
the cases that will be considered next, we assume Ai = 0.001, 
A 2 = .002, Pi = 53 dBm, P 2 = 33 dBm, B x = 1, B 2 = 1, 
^1 and ^2 are both exponential random variables with mean 
1, £1 = 3.8, £2 = 3.5, /3i = /?2 = dB, unless specified other 
wise. 

In our first example, we dwell in detail into our charac- 
terization of the SINR distribution as well as the coverage 
probability for arbitrary fading distribution that are differ- 
ent for different tiers. Figure [T] shows the plot of coverage 
probability versus the SINR threshold for various choices 
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Figure 2. 2-Tier network, Coverage Probability vs SINR threshold for various 
combinations of path-loss exponents 
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Figure 3. 2-Tier network, Coverage probability vs Tier 2 Bias Factor 

of fading distributions for each of the tiers. The first two 
curves (in the legend of Figure Q3 show the two-tier network 
coverage probability when the channel gains at both the tiers 
are exponential random variables with mean 1. In the next 
two curves in Figure Q] we depict two scenarios where the 
channel gains are log-normal random variables with zero mean 
and standard deviations a\ — 3.5 dB for the first tier and 
(72 = 4.65 dB for the second tier. The characterization of 
the coverage probability if fa ^ fa and if fa, fa dt dB 
were not known until now, for any chosen distribution of 
the channel gains, not even the exponential distribution. The 
results in this paper, in particular, Theorem [T] gives the 
coverage probability for all values of fa and fa. The last 
two curves in Figure Q] consider the channel gain to have 
exponential distributions with means /ii = 46.5 dB for the 
first tier and \i2 = 50.3 dB for the second tier. Notice that 
their coverage probability curves match exactly with those for 
the log-normal distributions. This is because, as illustrated in 



moment of the random fading factor of the k th tier, i.e. . 
If is a log-normal random variable with standard deviation 

' 2cr 2 \ 

and if is an exponential 



crfc, then E^fc = exp 1 
random variable with mean fik, the E^,' 



In this example, we have chosen ax, 02, /txi, and [12 in such 

a way that M i &i h 's are the same for k = 1, 2. 

In the next example, we show that we are now able to study 
the heterogeneous network for different path-loss exponents at 
different tiers. Notice that there is an improvement when the 
path-loss exponents are large. This is because the signal power 
decays faster with the distance, thereby causing lesser intercell 
interference. Notice from Corollary [T]that, when the path-loss 
exponents and the SINR thresholds are identical across the 
tiers, then the coverage probability is the same as that of a 
single-tier network with the same path-loss exponent and the 
SINR threshold, and further, coverage probability varies log- 
linearly with the SINR threshold (> dB). In Figure |2] notice 
that even when the path-loss exponents are not identical across 
the tiers, the coverage probability still has a log-linear behavior 
for SINR thresholds greater than dB. 

In the last example, we study the effect of varying the BS 
density (A2) and the bias factor (B2) of the second tier on 
the coverage probability, where the second tier has a greater 
path-loss exponent, to mimick a typical indoor environment 
situation. While Corollary Q] shows that the coverage proba- 
bility does not depend on the BS density of the tiers as well 
as the bias factors of the tiers when all the tiers have the same 
path-loss exponent and the SINR thresholds, when the path- 
loss exponents are different, Figure [3] shows that the coverage 
probability actually improves as we increase the bias factor 
of the second tier, and further it increases at a faster rate as 
the density of the second tier is increased. At the limits of 
the tier two bias factor, i.e. B2 — > 0, and B2 — > 00, the 
two-tier network essentially collapses to a single-tier network 
consisting of only the macrocell network, and the femtocell 
network, respectively. As mentioned previously, the single tier 
network is invariant to changes in the BS density and the bias 
factor, and as a result the curves are straight lines at these 
limits. 

The mathematical tools developed in this paper to study 
the heterogeneous network coverage probability are sufficient 
to characterize the average ergodic rate achieved at the MS, 
throughput and the per-tier traffic load, which are other im- 
portant metrics of interest for understanding the heterogeneous 
network, and these will be considered in detail elsewhere. 

VI. Conclusions 

In this paper, we study a heterogeneous cellular network 
consisting of K tiers, where each tier has its own BS density, 
BS transmission power and bias factor, path-loss exponent and 
channel gain with an arbitrary distribution, that is different for 
different tiers. For such a general model for the heterogeneous 
network, we develop mathematical tools based on stochastic 
geometry to characterize the distribution of the downlink SINR 
and the coverage probability at any given MS, where the MS 
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associates itself with the BS that has the SINR at the MS. 
Moreover, we have achieved a complete characterization of the 
SINR distribution and the coverage probability for all values 
of SINR thresholds, which has not been done before. 



network, corresponds to the BS that is nearest to the origin 
(MS) in the equivalent single tier network. As a result, SINR 
may be written in terms of the i?'s indexed in the ascending 
order, and we get © . 
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Appendix 

A. Proof for the Stochastic Equivalence Lemma 



B. Tail Probability of Desired BS Distance from MS 

From the proof of Lemma Q] the arrangement of the BSs of 
the fc th tier in the equivalent single-tier network is according 
to a non-homogeneous Poisson point process with BS density 

1 1—1 



networlcs^gHf" 



A fe f (P k B k ) 



> 0, 



Given a BS belonging to the k 



th 



Rk from the origin, then, R 



k = 



tier is at a distance 

(PkB^k)- 1 Rl* rep- 



resents the distance of the BS from the origin where the 
BS arrangement is according to a non-homogeneous 1-D 
Poisson point process with BS density function \( k > (r) = 



\ k f k {P k B k 
for each k 



t , r > 0, as long as < 00, 

1, 2, • • • , K. This is a consequence of the 
Mapping theorem [10 Page 18] and the Marking Theorem 
ifTUl Page 55] of the Poisson processes. Further, since the BS 
arrangements in the different tiers were originally independent 
of each other, the set of all the BSs in the equivalent 1-D 
non-homogeneous Poisson process is merely the union of all 
R's k, V k = 1, 2, • • • , K. By the Superposition Theorem 

iTTOl Page 16] of Poisson process, R (notice that it is not 
conditioned on k) corresponds to the distance from origin of 
BS arrangement according to non-homogeneous Poisson point 
process with density function A (r) = J\=i ( r ) > r — 0- 
In summary, we have converted the BS arrangement on a 
2-D plane of a if— tier network to a BS arrangement of the 
equivalent single-tier network along 1-D (positive x-axis), and 
further the SINR distributions of both these networks are also 
equivalent. Further, by our construction, the BS corresponding 
to the strongest received power at the MS, in the if —tier 



and Rf represents the distance of the nearest BS of this 
random process. Using the properties of the Poisson point 
process, P ({^ fe) > rj) = P ({N^ ([0,r]) = 0}) , where 

N^ ([0,r]) represents the average number of BSs in the 
interval [0, r] placed according to the Poisson point process 
with density A^ (r) , and this is equal to © • 

C. Desired BS Tier Probability Lemma 

The following sequence of equations provides the proof. 



P({/ = fe}) 



(a) 




< R 



{it® > R 



P»}) 



a ! tt(p /j b ; )~e*; 



dr. 



where (a) is obtained by noting that the desired BS belonging 
to the fc th tier is closer to the origin than the nearest BS (to 
origin) of the rest of the tiers, (b) is due to the fact that 

are independent random variables, and E R (k) [•] 
represents the expectation with respect to the random variable 
Rl , (c) is obtained by using Theorem [2}Equation Q to 
compute the probability of the event in (b) and to obtain the 

r (k) 

probability density function (p.d.f.) of R\ ' to evaluate the 
expectation, and finally, ([JJ is obtained by simplifying (c) . 

When {£/c}fLi = £, the integral in f4| simplifies to the 
form f t _ Q e~ at dt = i, which is rewritten in §5§ . 

D. Proof for the p.d.f. of the Desired BS Distance given Tier 
We first evaluate the probability of the event 
l^i > r, i = fcj . The steps are given in ([8]) , where 
(a) is obtained by noting that the desired BS belonging to 
the fc th tier is closer to the origin than the nearest BS (to 
origin) of the rest of the tiers, (b) rewrites joint probability 
in (a) in terms of the product of the marginal and the 
conditional probabilities, (c) is obtained by noting that the 
R<1\ V 1 = 1, 2,-- - , K are independent random variables, 
{d) is obtained by substituting for the tail probability events 
in (c) using Theorem [2} Equation ([3]) , and finally, (e) is 
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obtained by first evaluating the p.d.f. of R\ 
expressing the expectation as an integration. 
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As a result, 
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using ([3]) and 
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and finally, ([5]) is obtained by noting that /^j/Mfc) 



I = k 
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£. Proof for the Coverage Probability Theorem 

The sequence of equations in ^ , where (a) is obtained by 
using Lemma [1] and basic conditional probability properties, 
(6) is obtained by using [9, Theorem 1], using the BS density 
function specified in Lemma Q] for A (r) and 1 for the path- 
loss exponent, and (c) is obtained by exchanging the order 
of integrations in (6) . which is valid since the integrals are 
convergent. Upon simplifying, we get 



/ K JL 

i (wRi) = exp I -V :A' S 7T [PiBiR^j 



sl = l 



x ( 1 - a Fi ( -— ; l-—;iw 

Further, by evaluating the expectation in (c) by using Lemma 
|4] - Equation (J6j> , and simplifying, we get ([7]) . 

F. Simulation Method 

The fc th tier of the heterogeneous network with K tiers 
is identified by the following set of system parameters: 
(Afe, Pf., Bk, ££, /3fc) , where the symbols have all been 
defined in Section HH and k — 1, 2, • • • , K, where K is 
the total number of tiers. Now we illustrate the steps for 



simulating the heterogeneous network in order to obtain the 
SINR distribution and the coverage probability in the open- 
access cell association scheme. Assuming the MS to be at 
the origin, a single trial of heterogeneous cellular network 
arrangement in a cellular area with Rb as the boundary radius 
involves: 

1) Generating the random numbers iVfc ~ 
Poisson (Afc7ri?^) , which is the number of BSs of the 
k th that will be deployed in the trial. 

2) Generating Nk random variables according to a uniform 
distribution in the circular region of area ttR b , which rep- 
resents the location of the k th tier BSs corresponding to the 
trial. 

3) Computing the received power at the MS at the origin 
and computing the SINR as the ratio of the maximum of the 
received powers to the difference of the sum of all the received 
powers and the maximum received power. 

4) Record the index / which corresponds to the tier to which 
the desired BS belongs, for the tier. 

Repeat the same procedure T times. Typically, T is at 
least 50000. After this, we have an array containing the 
instantaneous SINRs and the tiers to which the desired 
BSs belonged, corresponding to the T trials. The tail prob- 
ability of SINR at a certain point, say r), is given by 

{# of trials where SINR > v \ , , , , ., 
7p '-, and the coverage probabil- 
ity of the MS in the heterogeneous network is given by 
{# of trials where i=k and siNR>p k \ 
- j, . 



